Magnetic Inelastic Dark Matter: Directional Signals Without a Directional Detector 
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The magnetic inelastic dark matter (MiDM) model, in which dark matter inelastically scat- 
ters off nuclei through a magnetic dipole interaction, has previously been shown to reconcile the 
DAMA/LIBRA annual modulation signal with null results from other experiments. In this work, 
we explore the unique directional detection signature of MiDM. After the dark matter scatters into 
its excited state, it decays with a lifetime of order 1 fis and emits a photon with energy ~100 keV. 
Both the nuclear recoil and the corresponding emitted photon can be detected by studying delayed 
coincidence events. The recoil track and velocity of the excited state can be reconstructed from 
the nuclear interaction vertex and the photon event vertex. The angular distribution of the WIMP 
recoil tracks is sharply peaked and modulates daily. It is therefore possible to observe the direc- 
tional modulation of WIMP-nucleon scattering without a large- volume gaseous directional detection 
experiment. Furthermore, current experiments such as XENON 100 can immediately measure this 
directional modulation and constrain the MiDM parameter space with an exposure of a few thousand 
kg ■ day. 
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I. INTRODUCTION 



Despite decades of direct detection efforts [T[ , the na- 
ture of dark matter interactions with regular matter re- 
mains elusive. The results from the DAMA/Nal and 
DAMA/LIBRA collaborations suggest that such inter- 
actions may be more intricate than originally expected. 
DAMA has observed an annual modulation in Nal crys- 
tals for the past decade 0-Q , with the expected phase for 
WIMP-nuclei interactions. There is no experimental evi- 
dence corroborating this signal. By now, it appears that 
the signal is not conventional spin-independent clastic 
scattering of WIMPs on nuclei. 

Among the quantitative explanations of DAMA, it is 
possible to take a few approaches. One method is to ex- 
ploit detector effects, such as channeling @-0|- Another 
is to introduce a dark matter model with more ingredi- 
ents (for example, [143 )• A possibility is that the dark 
matter preferentially scatters off the Nal used in DAMA, 
as opposed to the nuclei used in other direct detection 
experiments. In particular, we focus on the fact that io- 
dine is special in having both a relatively large mass and 
a relatively large magnetic moment [l4| . 

If dark matter has (weak) electromagnetic moments 
[Pol [lr| . it can interact through the charge and mag- 
netic dipole moment of the nuclei. For a summary of 
the interaction strengths for various nuclei used in direct 
detection experiments, see 17|. This type of interaction 



has been used to explain some recent direct detection 
results [13413' including the positive claim of DAMA. 
However, there are strong constraints from CDMS (23[ 
and XENON [H [H[ on this explanation of DAMA. 

Inelastic dark matter (iDM) takes advantage of the 
large iodine mass (26|. In iDM, there is an excited state 
of dark matter with mass splitting 5. It is assumed 
that dark matter couplings with nuclei are primarily off- 



diagonal, so that the WIMP is scattered into its excited 
state. This interaction only occurs if the dark matter has 
sufficient initial velocity. The minimum velocity v m i n for 
a WIMP to scatter with a nuclear recoil of energy En is: 
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where is the nucleus mass and m x is the WIMP 
mass. For splittings of 5 ~ 100 keV, experiments are only 
sensitive to the tail of the WIMP velocity distribution, 
leading to a much larger annual modulation than in the 
elastic case. Furthermore, scattering on heavier nuclei 
like iodine is preferred if m x is of order 100 GeV. 

The basic iDM model is now ti ghtl y constrained 12711 
by the latest results from CRESST [H, ZEPLIN-III [2g, 
XENON [Hi, and CDMS. It is possible to combine in- 
elastic scattering with yet more ingredients. For example, 
spin-dependent inelastic scattering is discussed in 



A. Magnetic Inelastic Dark Matter 



t.mpK 



We focus on magnetic inelastic dark matter (MiDM), 
because it has a unique and interesting directional sig- 
nature. Chang et al. [HI showed MiDM could explain 
both DAMA and other null results. The model takes ad- 
vantage of both the magnetic moment and large mass of 
iodine. In MiDM, the dark matter couples off-diagonally 
to the photon: 



C D (^) x*^uF^x + cc. 



(2) 



where the mass of x an d x* are s P u t by 5 ~ 100 keV. 
The off-diagonal coupling is natural if the dark matter 
is a Majorana fermion. The excited state has a lifetime 
t = ir/(5 3 fxi) ~ 1 — 10/iS, and emits a photon when it 
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TABLE I: In the first three (starred) rows, we give the best fit benchmark models of MiDM, with v eS c = 550 km/s and wo = 220 
km/s [lij . We also list parameters within the 90% CL region of each best fit value, for which the lifetime, r, can be a factor 
of a few larger or smaller. A is the average recoil track length. 77.15 is an estimate of the efficiency of XENON100 to detect 
delayed coincidence events, as described in Section III Bl The 'angular' rate is the rate for delayed coincidence events with a 
nuclear recoil in the energy range 10 — 80 keVr, followed by a photon with 5 keVee. This is obtained from multiplying the total 
rate by 77.15. We also show the expected number of nuclear recoil events for the published XENON100 non-blind analysis. 



decays. This short lifetime makes it possible to observe 
both the nuclear recoil and the emitted photon with a 
meter-scale detector. The two interaction vertices allow 
reconstruction of the excited state track. Both the ve- 
locity and angle can be measured, enabling directional 
detection even without a directional detector. 

A dark matter particle with a permanent electromag- 
netic dipole moment generally can be constrained by, 
e.g., gamma-ray measurements, the CMB, or precision 
Standard Model tests [3ll-[33| . However, the strongest 
bounds tend to come from direct detection experiments 
themselves, at least in the 100 GeV mass range. Fur- 
thermore, in MiDM, the inelastic nature of the interac- 
tion suppresses interactions with photons and baryons at 
low energies. If the dark matter is a composite particle, 
a low compositeness scale can also suppress annihilation 
to photons. 

There are some variants of the MiDM idea. In [34j . 
the parameter values were taken to be m x ~ 1 GeV and 
S ~ 3keV. The DAMA signal is produced by the emit- 
ted photon. This explanation evades constraints from 
other direct detection experiments because such low- 
energy electromagnetic events are typically rejected or 
not seen by other detectors. 

It is also possible that the dark matter couples to a new 
'dark' U(l), with gauge boson mass tua ^ p2l. l35l [3r| . 
Here the dark matter has a large dark dipole. If the dark 
gauge boson couples to regular electromagnetic currents, 
a sizable interaction with nuclei can be generated. How- 
ever, the decay rate of the excited state is suppressed 
because there is no direct interaction with the photon. 
While these ideas are interesting explanations of the 
DAMA signals, we do not consider them further because 
the excited state has a long lifetime. 

We study MiDM benchmarks, given in Table U which 
are good fits to the DAMA annual modulation signal 
flU ]. MiDM models with m x greater than ~ 300 GeV 
are severely constrained by ZEPLIN-III (29|, KIMS [37j], 
and XENON100 [U. 



The benchmarks are subject to form factor and veloc- 
ity distribution uncertainties {38U421 ] , especially for larger 
masses. The directional signal prediction can change 
wildly depending on the lifetime and rate. 

In order to explore the parameter space, we also con- 
sidered two extreme points within the DAMA 90% con- 
fidence level region found by [l4|, for each of the three 
masses. For the point with highest 6 and fi x , the ex- 
pectation for directional detection is better. The point 
with lowest S and /j x , which would not result in many 
delayed coincidence events, is in any case already tightly 
constrained by the XENON100 non-blind analysis. 

In this paper, we show that the current generation 
of direct detection experiments can observe a direc- 
tional signal from MiDM. For concreteness we focus 
on a XENONlOO-like detector, for two reasons. First, 
XENON100 will soon place strong constraints on the 
MiDM parameter space, making it the most relevant ex- 
periment to consider. Second, we wish to emphasize the 
feasibility of detecting a directional signal with experi- 
ments that are currently running. 

We compute the distribution of recoil track angles 
and velocities from MiDM benchmarks. The sensitivity 
of XENON100 to the MiDM parameter space depends 
strongly on the lifetime of the excited state. For the 
benchmark lifetimes of ~ 1 — 10/js, XENON100 can mea- 
sure the directional modulation at high significance and 
obtain sharp constraints on the parameter space with 
just tens of events. This is achievable with around 5000 
kg • day in the energy range 10 — 80 keVr. 



II. DIRECTIONAL DETECTION 

Directional detection can clearly test whether any sig- 
nal comes from WIMP interactions [43j]. Due to the 
Earth's motion in the Galaxy, there is a "WIMP wind" 
which is opposite the motion of the Earth. There is a 
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daily modulation in the angle of recoil tracks in the lab 
frame. This modulation depends only on the rotation of 
the Earth relative to the WIMP wind, and can be disen- 
tangled from the daily rotation of the Earth with respect 
to the Sun. The experimental directional detection ef- 
fort focuses on measuring the nuclear recoil track with 
large- volume, gaseous detectors (43 - l46j . 

Angular information is a particula rly pow erful discrim- 
inant of WIMP scattering for iDM [13,111] ■ Because in- 
clastic interactions have a high velocity threshold, the 
angular distribution of the nuclear recoil tracks is sharply 
peaked in the direction of the WIMP wind. There is a 
kinematic constraint on the recoil angle of the nucleus: 



{COS "f) max (Er) 



ve 



(3) 



Here 7 is the angle between the velocity of the Earth and 
the recoil velocity in the Earth frame, ve is the Earth's 
velocity in the Galactic frame, and v esc is the Galactic 
escape velocity from the Solar neighborhood. For typical 
iDM models considered in the literature, 7 is constrained 
to be within ~100 degrees of the WIMP wind [47|. How- 
ever, because the signal goes to zero at the bound in 
Eq. [3l the precise location of this kinematic constraint 
can be difficult to pinpoint. 

MiDM has better directional detection prospects at 
XENON100, compared to directional detection of iDM. 
Current directional detectors focus on spin-dependent 
scattering and use light targets such as CF 4 [49H52) . 
Thus, they would not see inelastic scattering events. In 
the MiDM case, there is also much more event informa- 
tion and thus more sensitivity to the parameter space. 
One can measure both the velocity («*) and the angle 
(cos 7*) of the WIMP recoil track. Once again, this re- 
coil angle is with respect to the Earth's motion. The 
tracks are sharply peaked in angle opposite the motion 
of the Earth. 

For the WIMP recoil angle, there is also an energy- 
dependent maximum recoil angle, which we give in 
Sec. Mil The most important bound is on the WIMP 
recoil velocity, 



(E R (m N /m x - 1) - 5) 



y/2m N E R 



(4) 



Here the signal peaks near the kinematic bound because 
most events occur near the threshold velocity in Eq. Q] 
Thus having information on both and E R is an ex- 
tremely sensitive probe of the model parameters. There 
is a remaining degeneracy: if S and m x are shifted in 
opposite directions, the bound can remain roughly the 
same. However, one can fit S separately from the spec- 
trum of the nuclear recoils, and from the energy of the 
emitted photons. 

There is also a maximum velocity for the excited state, 



vr x {E R ) = ^/(v E + v esc y-2(E R + 5)/m x , (5) 
but the rate is exponentially suppressed at this bound. 



A. XENON100 

We model directional detection in XENON100 with a 
simplified XENONlOO-like experiment. XENON10 0, 
HH had 316.4 kg- day of data in the energy range 4.5- 
75 keVr. XENON 100 has a 40 kg fiducial mass, at even 
lower backgrounds. The initial 170 kg- day non-blind run 
already constrains the MiDM parameter space (at low 6) . 

The XENON100 detector is a cylinder, with a radius of 
15.3 cm and a height of 30.6 cm. The fiducial volume has 
a radius of 13.5 cm and height of 24.3 cm. The primary 
scintillation (SI) and ionization (S2) signals of an event 
are measured. For more details, see j54|. The S2 signal 
is observed 15-140 /is after the SI signal, for events in 
the fiducial volume. 

The signature of MiDM is two SI signals separated by 
roughly .5 fis in time, followed at least 15 (is later by two 
S2 signals. The photon event is identified from the second 
SI signal and an S2 peak with energy of ^100 keVee. 
At 100 keVee, a photon is clearly distinguishable from 
a nuclear recoil by S2/S1. The other event should be 
consistent with a nuclear recoil. The time separation of 
the two S2 signals depends on how the WIMP recoils 
along the cylinder axis, z. In XENON10, events with 
multiple S2 events at different z positions were rejected. 

We refer to the track connecting the two events as the 
decay track. Events can be localized to a 3D spatial 
resolution of 3 mm (though the absorption length for the 
100 keV photon may blur this) and timing resolution of 
10 ns. Meanwhile, the track should be at least 10 cm long. 
This makes it possible to measure direction and velocity 
of the decay track to an extremely high accuracy. The 
head-tail discrimination of the track can be determined 
using timing information and the S1/S2 ratio. 

We wish to obtain the \* recoil track from the decay 
track. However, because the photon can travel up to ~1 
cm after emission, this introduces systematic uncertain- 
ties. The observed decay track can be blurred by a few 
degrees, relative to the x* recoil track direction. This 
also introduces an uncertainty in the velocity of the x* 
of roughly 10%. 

There are some specific event geometries that can re- 
sult in more ambiguous events. For example, it could be 
difficult to resolve the two S2 signals if the decay track 
is perpendicular to the z axis. Then the two S2 signals 
arrive at nearly the same time. S2 signals generally have 
a time width of ~ l/is and the PMT spatial resolution 
is only ^2.5 cm. However, because the drift velocity is 
2mm//is, this is a small fraction of the total solid angle. 

Thus directional events are in principle detectable at 
XENON100. The background for such delayed coinci- 
dence events with both a nuclear recoil and a photon 
of energy ~ 100 keV should be extremely low. There 
are other 'mixed' delayed coincidence events from Bi and 
Kr contamination, and excitation of metastable states of 
Xe [HI]. However, these have very different energies and 
decay times. It may be possible to extend the fiducial 
volume when searching for directional events. 
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FIG. 1: The efficiency r\ for the best-fit benchmarks from 
TableU Ro is the size of a spherical detector. We approximate 
the XENON100 fiducial volume as a sphere with radius Rq = 
0.15 m, marked by the vertical black line. The thinner lines 
show the corresponding results with highest S, within the 90% 
CL region of the best fit. 



B. Detector Efficiency 



L = v*t' . Here we neglect the smearing arising from the 
mean free path of the emitted photon, since the emission 
is isotropic. There is also no dependence on t or recoil 
angle in this approximation. Then the expression for 
efficiency above can be simplified to 

dLg(L) J ^_^H(R,L) (6) 

where L is the recoil length. The recoil length distribu- 
tion g{L) is 

g{L)= ( dVr iM exp (_±\ (7) 

J «*T V V * T J 

where now /(i>*) is the distribution for v*, not u*. A good 
approximation is g{L) = exp(— L/X)/X, where A = (v*)t 
is the average recoil length. Typical A values are given 
in Table [J 

We approximate the XENON100 detector as a sphere. 
The fiducial volume has radius Rq = 15 cm, with effi- 
ciency 77.15. Results are shown in Fig. [1] The precise 
efficiency depends on specifics of the detector, and must 
take into account the effects mentioned in Sec. Ill Al 



The typical decay length is 1 — 10 m in these mod- 
els, relatively large compared to XENON100. Thus the 
WIMP can recoil inside the detector volume, but decay 
outside the detector Jg3|. The effective exposure for de- 
layed coincidence events is, in general, lower than the ex- 
posure for nuclear recoils because of this geometric effect. 
Here we compute the detector efficiency, as a function of 
typical detector size, for the MiDM benchmarks. 

The efficiency is 



d 3 



df 



-t'h 



fa 
V 



H(t,x,d,t') 



The term in brackets comprises detector effects. The spa- 
tial integral is over the detector volume. H(t,x,d,t') is 
the efficiency for observing a WIMP decay, given that a 
nuclear recoil was observed. This depends on the time 
of the year t, the location of the WIMP-nucleus inter- 
action inside the detector, x, the decay vector, d, and 
the WIMP decay time (coincidence time), t'. Whether 
a given WIMP decay track is located inside the detector 
depends on the orientation of the detector with respect to 
the Earth's velocity, the decay vector, and the efficiency 
for the particular event geometry. 

The astrophysics and particle physics is captured by 
the integral over t' and v„ . r is the lifetime of the excited 
state. The distribution of recoils depends on the WIMP 
recoil velocity distribution, /(??*), and the decay time 
distribution. We assume that iT* is defined with respect to 
the Earth's velocity vector so that f(v*) does not depend 
on detector orientation. 

For the calculation below, we model the detector as 
a single sphere of size i?o- We assume that H depends 
only on the interaction position x and the decay length 



III. RECOIL SPECTRUM 

There are two electromagnetic scattering channels for 
magnetic dark matter: dipole-dipolc and dipole-chargc. 
In the dipolc-dipole scattering case, the dark matter in- 
teracts with the magnetic moment of the nucleus. The 
matrix element is 

\M\ 2 2 [ ^nuc\ 2 [ V X \ 2 Sn+± 
— — = K)7T(X TTLjsJ ( — , (8 

62irmisrm x V e / V e / 6b n 

for each isotope. We sum over all isotopes, weighted by 
their abundances fvj\ . There is, in general, also a factor 
of (S x + 1)/(3S X ) for the spin of the dark matter. We 
take S x = 1/2. 

In the dipolc-charge scattering case, the dark matter 
interacts with the electric charge of the nucleus. The 
matrix clement is 



32irm^m 2 c 



En \ e ) \2mjv m xJ 



(9) 



where v is the initial velocity of the WIMP in the lab 
frame. We have again assumed S x = 1/2. 

The differential scattering rate for measuring both nu- 
clear recoil energy and WIMP recoil track is 



dR 



r]N T p x 



dEndv^dx* 



d 3 v f(v + ve) v 



da 



dEfidv*dx* 
(10) 

where we have abbreviated x* — cos 7*. The three- 
dimensional WIMP velocity distribution is given by f(v). 
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FIG. 2: Differential rates dR/(dE R dcos-y t dv t ) for (a) m x = 140 GeV MiDM benchmark and (b) m x = 70 GeV MiDM 



benchmark. The 77K 



300 GeV benchmark looks similar to the rrn 



140 GeV benchmark. In each case we show three 



two-dimensional distributions, where we have integrated over the third variable. All rates are computed assuming scattering 
on Xe, and benchmarks are given in Table U The differential rate is normalized so that the total rate is unity. 



Nt is the number of target nuclei per kg and p x is 
the local WIMP energy density, which we fix to be 0.4 
GeV/cm 3 " 



55|. 



As in [47J , we expand da and change variables to v' = 
v + Ve- The trivial integral over v 1 imposes the condition 

v = q/m x + v* + ve- (11) 

q is the recoil momentum of the nucleus. The resulting 
differential rate is 



dR 



F 2 [E TI } 



5 d)l 9- 



2m, 



q-v*-E R -5). (12) 



F 2 [Eft\ is a nuclear form factor which depends on the 
type of interaction. 

For a xenon target, dipole-charge scattering, Eq. [9l 
dominates. For this we use the standard nuclear Helm 
form factor. Dipolc-dipole scattering, Eq. [51 is roughly 
20% of the total rate. To calculate dipolc-dipole scatter- 
ing a magnetic moment form factor is necessary. The nu- 
clear magnetic moment receives contributions from both 
spin and angular momentum. We use the spin form fac- 
tor from [561 ] . The angular momentum component is 
~ 20—30% at zero momentum for Xc. Since dipolc-dipole 
scattering is already subdominant for Xe, and since we 
do not have accurate angular momentum form factors, 



we approximate the entire magnetic moment form factor 
with the spin component. 

We now specialize to the case where f(v) is a normal- 
ized, truncated Maxwell-Boltzmann distribution, with 
v esc = 550 km/s [13] and vq = 220 km/s. We assume 
ve = 240 km/s on average and label the normalization 
factor of the distribution as n(i>o, v esc ). The result is 



dR 



dEftdv*dx* 



7]N T p x v* | -M| 2 
m x 327rmjvm 2 , 



F 2 [E R }Q(l~\x q \) 



m (is) 



and 



(14) 



n(v ,v esc ) 
with the following definitions: 

_ {E R (m N /m x - 1) - S) 

Xn 

qv* 

(v 1 ) 2 =v E + q 2 /m x + v 2 + 2v E v*x* + 2x q v*q/m x 
+ 2vEq/m x (x q x* + 1 — x 2 \J\ — x 2 cos < 

An upper bound on a;* can be extracted from setting 
v' = v esc , with cos0 = — 1. The bound depends on both 
V* and E R . 

Finally, the matrix elements are given in Eq. [8J or in 
Eq. O Note that in the dipole-charge scattering case we 
need to replace v in Eq. [9] using the energy conservation 
relation, m x v 2 = 2E R + m x v 2 + 25. 

The normalized total rate spectrum of several bench- 
mark models is shown in Fig. [5] 
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FIG. 3: Exposure to obtain a 5a measurement of (cos 7.) 90% of the time the experiment is conducted on Earth. The energy 
range of the experiment is 10 — 80 keVr. (IRbg is the background rate; a XENON maximum background rate is indicated by 
the solid vertical line. The bands shown give the exposures necessary as the rates modulate throughout a year. We show (a) 
three mass benchmarks and (b) three mass benchmarks in the 90% CL region with highest delta, from Table [I] 



IV. SENSITIVITY 



B. Parameter Estimation 



XENON100 is collecting several thousand kg- day of 
exposure. We assume a total exposure of 5000 kg- day 
on a 40.6 kg fiducial target, in a nuclear recoil energy 
range of 10-80 keVr. This is consistent with scaling up 
the results from XENON10 and with preliminary results 
reported by XENON100. 

For the best-fit parameters listed in Table [I] this would 
imply a minimum of ~ 100 nuclear recoils observable by 
XENON100. Only ~ 10 delayed coincidence events are 
expected, due to the small size of the detector relative 
to the average recoil track length. Despite these low effi- 
ciencies, a study of the delayed coincidence events is still 
vastly more informative in two ways: (a) it establishes a 
directional signal correlated with the WIMP wind, and 
(b) it is much more sensitive to the parameter space. 



A. Directional Detection 

We first determine the exposures required to establish 
a correlation with the WIMP wind. The average nuclear 
recoil angle with respect to the Earth's motion, (cos 7), is 
a robust model-independent statistic for directional de- 
tection [13, [H, H3|. Here we use (cos 7*), the WIMP 
recoil angle with respect to the Earth's motion. Because 
of the rotation of the Earth, on average (cos 7) or (cos 7*) 
should be consistent with for standard backgrounds. 

Because XENON100 has excellent spatial resolution, 
we assume that the recoil track angle can be determined 
to 10 degrees. We compute the exposures required to 
obtain a 5a result for (cos 7*) 90% of the time. We 
allow for a uniform (isotropic) background, though the 
XENON100 background should be negligibly low. The 
results are shown in Fig. [3J The required exposures 
roughly correspond to a minimum of 16 events at zero 
background. 



The predicted rate for delayed coincidence events at 
XENON100 is only a few counts per 1000 kg ■ day. How- 
ever, the additional recoil track information makes it pos- 
sible to obtain an excellent measurement of the model 
parameters. 



We perform a likelihood analysis, as described in |47j . 



over the parameter space of m X7 5, and fx x . We compute 
the (relative) log likelihoods for £ kg- day on Xe, with 
nuclear recoil energy range 10 — 80 keVr. We neglect the 
effects of imperfect angular and energy resolution. (The 
XENON100 energy resolution is ~ 10% in this energy 
range, and we estimate an angular resolution of 10 de- 
grees.) The log likelihood is 

laCtotip) = £ J dx^i(x;p Q )\nfj,(x;p) - fi(x',p)\ (15) 

where p refers to (m x ,5,fi x ) and x refers gcncrically to 
the event space of cither E R or (E R , u„, cos 7*). f>o are 
the true model parameters. /1 is the rate for parameters 
p. If there is only nuclear recoil energy information, 



h{Er\p) 



dR 
dE~R 



(E R ;p) + dR B a, 



(16) 



in cpd/kg/kcVr for parameters p. We assume the back- 
ground rate, dR R c is known and negligibly small. 
If there is directional information, 



n(E R ,v*,x*;p) = 77.15(f)) 



dR 



dE R dv*dx.. 



dR 



BG 



dv*dx* 



where x* = cos 7*. 77.15(f)) is the efficiency, for parameters 
p, at XENON100. 

In Fig. [4] we show the sensitivity to MiDM parameters 
if (1) only nuclear recoil information is used and (2) if 
only delayed coincidence events arc considered for 5000 
kg • day. We show confidence levels of (68, 90, 95, 99, 
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FIG. 4: Confidence levels for determining m x , 5, and fj, x , marginalized over the third parameter for each two-dimensional slice. 
We assume an exposure of 5000 kg- day on Xe in the energy range 10 — 80 keVr. The plots show sensitivity to the MiDM 
parameter space, using only delayed coincidence data, for the (a) m x — 70 GeV benchmark, and (b) m x = 140 GeV benchmark, 
(c) shows the sensitivity using only nuclear recoil events, for the m x = 140 GeV benchmark. The directional information is a 
better test of m x and fj, x . The case with m x — 300 GeV looks similar to m x = 140 GeV. 



and 99.9%). We neglect the Earth's velocity about the 
Sun since a livetime of order a year is needed for 5000 
kg • day. 

Despite the reduction by a factor of 10-50 in events, the 
directional data is a much stronger constraint on m x and 
fi x . S can also be determined from the Er data or the 
photon energies. In Fig. [5] we show the sensitivity to m x 
and \x x for the m x = 140 GeV benchmark, assuming that 
6 is already known. The directional information breaks 
the degeneracy in m x and [i x when only nuclear recoil 
information is used. 



C. Measurement of Both Recoils 



So far, we considered measurement of the WIMP recoil 
velocity vector from delayed coincidence events. With a 
gaseous directional detector, it is also possible to obtain 
the recoil angle of the nucleus. Then m x and 8 are highly 
constrained. For such events there are 4 equations and 
5 unknowns: m Xl 6, and v. However, one can obtain 6 
from the energy peak of the coincident photons. Then it 
is possible to measure the WIMP mass and velocity with 
just 1 WIMP scattering event. The mass is determined 
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FIG. 5: Confidence levels for determining m x and fi x using 
WIMP recoil tracks, assuming S is already measured from the 
photon energies or nuclear recoil spectrum. Here we take the 
m x = 140 GeV benchmark and assume an exposure of 5000 
kg- day on Xe in the energy range 10 — 80 keVr. 



V. CONCLUSIONS 

The magnetic inelastic dark matter model has an in- 
teresting and previously unstudied signature at direct de- 
tection experiments: a delayed coincident photon with 
energy 5. Observation of such photons would also allow 
current direct detection experiments to become excellent 
directional detectors. 

Motivated by the MiDM setup, we studied several 
benchmark model parameters that can fit the com- 
bined DAMA/Nal and DAM A/LIBRA data. Given the 
rapidly improving constraints from other experiments, 
we feel that MiDM is currently the best hope for a 
dark matter interpretation of DAMA - and it predicts a 
low-background signature detectable with current experi- 
ments. 

With 5000 kg- day of exposure, XENON100 can detect 
the angular modulation of the recoils and determine 
the MiDM model parameters. While we focused on 
benchmarks in MiDM, we emphasize that such a delayed 
coincidence signal is worth searching for in general. 
Such events, if found, carry much more information 
than simple nuclear recoils, and would provide more di- 
rect access to the WIMP velocity distribution in our halo. 



by the following equation: 

2m N E R 

m x = — — , — . (17) 

2(6 + En) — \/2mNER q ■ v* 

Since q and v* are measured, the initial WIMP velocity v 
is then fixed by momentum conservation. A direct mea- 
surement of the WIMP velocity distribution is then also 
possible. 
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